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Abstract 

This  study  presents  molecular  orbital/density  functional  theory  (MO/DFT)  calculations  of  the  electronic  structure,  vibra¬ 
tional  frequencies,  and  equilibrium  isotope  fractionation  factors  for  iron  desferrioxamine  B  (Fe-DFO-B)  complexes  in  aque¬ 
ous  solution.  In  general,  there  was  good  agreement  between  the  predicted  properties  of  Fe(III)-DFO-B  and  previously 
published  experimental  and  theoretical  results.  The  predicted  fractionation  factor  for  equilibrium  between  Fe(III)-DFO-B 
and  Fe(III)-catecholate  at  22  °C,  0.68  ±  0.25%o,  was  in  good  agreement  with  a  previously  measured  isotopic  difference 
between  bacterial  cells  and  solution  during  the  bacterial-mediated  dissolution  of  hornblende  [Brantley  S.  L.,  Liermann  L. 
and  Bullen  T.  D.  (2001)  Fractionation  of  Fe  isotopes  by  soil  microbes  and  organic  acids.  Geology  29,  535-538].  Conceptually, 
this  agreement  is  consistent  with  the  notion  that  Fe  is  first  removed  from  mineral  surfaces  via  complexation  with  small  organic 
acids  (e.g.,  oxalate),  subsequently  sequestered  by  DFO-B  in  solution,  and  ultimately  delivered  to  bacterial  cells  by  Fe(III)- 
DFO-B  complexes.  The  ability  of  DFO-B  to  discriminate  between  Fe(III)  and  Fe(II)/Al(III)  was  investigated  with  Natural 
Bond  Orbital  (NBO)  analysis  and  geometry  calculations  of  each  metal-DFO-B  complex.  The  results  indicated  that  higher 
affinity  for  Fe(III)  is  not  strictly  a  function  of  bond  length  but  also  the  degree  of  Fe-O  covalent  bonding. 

©  2008  Elsevier  Ltd.  All  rights  reserved. 


1.  INTRODUCTION 

Iron  is  an  essential  nutrient  for  the  majority  of  organ¬ 
isms  and  functions  in  a  variety  of  cell  reduction  and  enzy¬ 
matic  processes.  Despite  the  high  natural  abundance  of 
iron,  its  availability  to  plants  and  microorganisms  under 
aerobic,  neutral  pH  conditions  is  limited.  Iron  limitations 
primarily  arise  from  the  low  solubility  of  Fe(III),  which  at 
circumneutral  pH  is  ~10-18  M  with  respect  to  most 
Fe(III)-bearing  minerals  (Kiss  and  Farkas,  1998).  Because 
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of  the  difference  between  the  biological  demand  for  Fe 
and  its  availability,  accumulation  of  Fe  by  organisms  signif¬ 
icantly  impacts  its  cycling  in  natural  systems  and  Fe  can  be 
a  limiting  factor  for  biological  productivity  (Martin  and 
Fitzwater,  1988). 

Siderophores,  an  important  class  of  organic  acids  with 
large  complexation  constants  for  Fe,  are  produced  by  sev¬ 
eral  organisms  in  order  to  overcome  iron  deficiencies  (Wie- 
derhold  et  al.,  2006).  Due  to  their  exceptionally  high  affinity 
for  Fe,  siderophores  complex  Fe(III)  by  extracting  it  from 
Fe-bearing  minerals  or  aqueous  Fe(III)  complexes.  Fe-sid- 
erophore  complexes  can  then  be  delivered  to  an  organism 
for  uptake  through  a  variety  of  cellular  recognition  and 
transport  mechanisms  (Albrecht-Gary  and  Crumbliss, 
1998;  Boukhalfa  and  Crumbliss,  2002).  Subsequent  release 
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of  Fe,  at  or  through  an  organism’s  cell  membrane,  allows 
the  siderophore  to  scavenge  for  additional  Fe. 

Complexation  of  Fe(III)  by  naturally  synthesized  sidero- 
phores  is  efficient  because  of  the  remarkably  large  associa¬ 
tion  constants  for  Fe(III)-siderophore  complexes. 
Association  constants  for  Fe(III)-siderophore  complexes 
range  between  1023  to  1052  (Hernlem  et  al.,  1996;  Albr- 
echt-Gary  and  Crumbliss,  1998).  Chemical  specificity  for 
Fe(III)  is  readily  achieved  because  association  constants  be¬ 
tween  siderophores  and  other  di-  and  tri-valent  metals  are 
relatively  smaller.  For  example,  the  association  constants 
for  desferrioxamine  B  (DFO-B)  (a  natural  siderophore  syn¬ 
thesized  by  the  soil  actinomycete  Streptomyces  pilosus)  with 
different  metals  are  as  follows:  Fe(III)  =  1042’33,  Ga(III) 
=  1038-96,  Al(III)  =  1036  11,  Ni(II)  =  1027-66, 
Cu(II)  =  1023'98,  and  Zn(II)  =  1020'40  (Anderegg  et  al., 
1963;  Evers  et  al.,  1989;  Kurzak  et  al.,  1992;  Hernlem 
et  al.,  1996).  This  chemical  specificity  is  important  in  envi¬ 
ronments  such  as  soils,  where  Al(III)  concentrations  may  be 
comparable  to  Fe(III)  concentrations.  Because  Al(III)  is 
phytotoxic  (Kochian,  1995,  and  references  therein),  the 
ability  to  uptake  Fe(III)  and  exclude  Al(III)  is  critical  to 
plant  health. 

Cheah  et  al.  (2003)  reported  that  the  dissolution  of  goe- 
thite  in  the  presence  of  both  oxalate  and  DFO-B  was  en¬ 
hanced  relative  to  experiments  involving  just  one  of  these 
two  organic  acids;  the  authors  suggested  that  Fe  transport 
from  a  mineral  surface  to  a  bacterial  cell  involves  a  con¬ 
certed  multi-step  process.  First,  oxalate  extracts  Fe  from 
the  mineral  surface  and  Fe-oxalate  complexes  diffuse  into 
solution.  Second,  DFO-B  scavenges  Fe  from  Fe-oxalate 
complexes,  freeing  oxalate  to  return  to  the  mineral  surface. 
Finally,  Fe-DFO-B  complexes  transport  Fe  to  bacterial 
cells  and  release  of  Fe  allows  DFO-B  to  scavenge  additional 
Fe. 

Other  studies  have  measured  the  fractionation  of  Fe  iso¬ 
topes  during  leaching  of  Fe  from  hornblende  (Brantley 
et  al.,  2001)  and  goethite  (Brantley  et  al.,  2004).  The  former 
of  those  two  studies  reported  that  the  amount  of  Fe  isotope 
fractionation  between  hornblende  and  solution  depended 
on  the  Fe-ligand  affinity.  That  study  also  noted  an  isotopic 
fractionation  between  Fe  that  remained  in  solution  and  Fe 
that  was  incorporated  into  bacterial  cells  (0.85  and 
1.17  ±0.20%o).  In  these  isotope  studies,  a  major  challenge 
was  quantifying  and  discriminating  between  the  amount 
of  fractionation  caused  by  equilibrium  effects  in  solution 
and  by  kinetic  effects  at  mineral-solution  and  solution-cell 
interfaces.  The  first  measurements  of  fractionations  be¬ 
tween  organic-bound  and  inorganic  Fe  in  solution  recently 
appeared  in  the  literature  (Dideriksen  et  al.,  2008),  and  this 
study  aims  to  re-produce  the  fractionations  reported  in  that 
study. 

The  process  of  separating  ligand-bound  Fe  from  aque¬ 
ous  Fe  can  lead  to  isotopic  fractionation  (Wiederhold 
et  al.,  2006).  Thus,  validation  and  interpretation  of  experi¬ 
mental  findings  would  be  greatly  assisted  by  the  develop¬ 
ment  of  models  that  can  accurately  predict  isotope 
fractionation  factors.  Isotope  fractionation  factors  have  al¬ 
ready  been  predicted  and  measured  for  Fe(II)-(H20)6  and 
Fe(III)-(H20)6  (Johnson  et  al.,  2002;  Welch  et  al.,  2003; 


Jarzecki  et  al.,  2004;  Anbar  et  al.,  2005;  Domagal-Goldman 
and  Kubicki,  2008),  and  have  been  predicted  for  Fe(II)- 
oxalate  and  Fe(III)-oxalate  (Domagal-Goldman  and  Kub¬ 
icki,  2008).  However,  fractionation  factors  for  siderophore 
complexes  have  not  been  previously  predicted. 

The  objectives  of  this  study  were  to  examine  the  high  affin¬ 
ity  and  specificity  of  DFO-B  for  Fe(III)  and  to  predict  the 
equilibrium  Fe  isotope  fractionations  associated  with  the 
Fe-DFO-B  complex  in  solution.  First,  we  present  the  effects 
of  different  methodologies  on  the  predicted  structure  of 
DFO-B  complexes  because  Fe-O  bond  distances  influence 
both  Fe  isotope  fractionation  factors  and  Fe  affinities.  We 
then  examine  the  chemical  specificity  for  Fe  using  NBO  (Nat¬ 
ural  Bond  Orbital)  analysis  (Weinhold  and  Landis,  2001)  of 
Fe(III)-DFO-B,  Fe(II)-DFO-B,  and  Al(III)-DFO-B,  where¬ 
by  the  electron  occupancies  of  Fe-O/Al-O  bonds  and  the 
partial  charges  of  the  Fe/Al  atoms  are  predicted.  We  also 
present  predicted  vibrational  frequencies  for  Fe(III)-DFO- 
B  and  Fe(II)-DFO-B  complexes,  since  these  predictions  af¬ 
fect  the  accuracy  of  calculated  equilibrium  isotope  constants 
for  Fe(III)-DFO-B  and  Fe(II)-DFO-B. 

2.  METHODS 
2.1.  Complexes  studied 

The  siderophore  examined  in  this  study  was  DFO-B,  the 
most  extensively  studied  siderophore  with  respect  to  min¬ 
eral  dissolution  of  the  approximately  500  known  sidero¬ 
phore  structures  (Wiederhold  et  al.,  2006).  DFO-B  is  a 
linear  trihydroxamic  acid  composed  of  1,5-diaminopentane 
and  succinic  acid  residues.  DFO-B  has  four  acidic  protons 
(three  hydroxamic  acids  and  one  terminal  amine),  is  posi¬ 
tively  charged,  and  is  soluble  in  water  and  alcohols  (Borgias 
et  al.,  1989).  DFO-B  coordinates  Fe(III)  through  six  oxygen 
donors,  forming  three  chelate  rings  in  a  distorted  octahe¬ 
dral  symmetry.  The  resulting  Fe(III)-DFO-B  hexadentate 
complex  has  16  possible  geometrical  and  optical  isomers 
that  have  small  energy  differences  (Leong  and  Raymond, 
1975;  Borgias  et  al.,  1989;  Dhungana  et  al.,  2001). 

Previous  investigations  of  hydroxamate  complexes  at  the 
goethite-water  interface  have  been  conducted.  Holmen  and 
coworkers  (Holmen  and  Casey,  1996;  Holmen  et  al.,  1997) 
used  acetohydroxamic  acid  (aHa)  as  an  analog  for  high 
molecular  weight  hydroxamate  siderophores  and  per¬ 
formed  Hartree-Fock  calculations  to  investigate  various 
conformations  of  aHa.  Edwards  et  al.  (2005)  used  the 
B3LYP/6-311G(d)  method  to  predict  the  vibrational  fre¬ 
quencies  of  Fe-aHa  complexes.  To  date,  isotope  fraction¬ 
ation  factors  have  not  been  predicted  for  Fe-DFO-B. 
Under  circumneutral  pH  conditions,  Fe(III)-DFO-B  has  a 
+  1  charge  arising  from  the  terminal  amine  group.  However, 
in  the  crystal  structure,  the  +1  charge  is  balanced  by  a 
counter-ion.  For  the  gas-phase  calculations,  we  used  a  neu¬ 
trally  charged  Fe(III)-DFO-B  complex  in  order  to  minimize 
intra-molecular  H-bonding  interactions  that  do  not  exist  in 
solution  (Section  3.2).  For  the  solution-phase  calculations, 
we  used  the  protonated  Fe(III)-DFO-B  complex  because 
inclusion  of  solvation  prevented  the  formation  of  artificial 
intra-molecular  H-bonds  (Fig.  1). 
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2.2.  Quantum  chemical  calculations 

Molecular  orbital/density  functional  theory  (MO/DFT) 
calculations  were  performed  with  Gaussian  03  (Frisch  et  al., 
2004).  Specifically,  we  employed  the  unrestricted  B3LYP 
method,  comprised  of  Becke’s  (Stephens  et  al.,  1994;  Becke, 
1997)  3-parameter  nonlocal-exchange  functional  (B3),  and 
Lee,  Yang,  and  Parr’s  (Lee  et  al.,  1988)  gradient-corrected 
correlation  functional  (LYP).  The  B3LYP  method  has  been 
shown  to  accurately  predict  the  molecular  structures  and 
vibrational  frequencies  of  transition  metal  complexes 
(Koch  and  Holthausen,  2001),  including  Fe-aHa  and  an 
analog  for  DFO-B  (Edwards  et  al.,  2005). 

The  local  bonding  environment  of  Fe(III)-DFO-B,  as 
predicted  by  the  MO/DFT  calculations,  was  compared  to 
experimental  single-crystal  X-ray  diffraction  data  (Dhung- 
ana  et  al.,  2001)  and  EXAFS  measurements  of  the  solvated 
complex  (Edwards  and  Myneni,  2005).  MO/DFT  predic¬ 
tions  are  presented  with  the  following  basis  sets:  6-3 1G, 


3 

6-31G(d),  6-31+G(d),  6-31  lG(d),  6-311G(d,p),  and  6- 
31 1+G(d).  The  geometries  were  predicted  from  energy  min¬ 
imizations  performed  without  symmetry  constraints  that 
could  have  potentially  biased  the  predictions.  Frequency 
calculations  were  subsequently  performed  to  verify  that  a 
minimum-energy  structure  was  located  (i.e.,  no  imaginary 
frequencies). 

IR  frequencies  are  reported  for  gas-phase  Fe(III)-DFO- 
B  and  Fe(II)-DFO-B  complexes,  using  the  B3LYP/6- 
31  lG(d,p)  method.  A  similar  method  (B3LYP/6-311G(d)) 
accurately  predicted  the  IR  frequencies  for  Fe(III)-triscate- 
chol  and  Fe(III)-trisoxalate,  two  complexes  with  octahe¬ 
dral  coordination  similar  to  Fe(III)-DFO-B  (Domagal- 
Goldman  and  Kubicki,  2008).  Frequencies  were  scaled  by 
a  factor  of  0.956  in  order  to  account  for  errors  associated 
with  the  use  of  a  harmonic  oscillator  approximation  (Hehre 
et  al.,  1985)  as  well  as  the  inaccuracies  associated  with  the 
DFT  method.  This  scale  factor  was  developed  via  modeling 
of  smaller  Fe-organic  complexes  (Edwards  et  al.,  2005),  has 


d 


Fig.  1.  Optimized  geometries  for  protonated  gas-phase  Fe-DFO-B  (a);  de-protonated  gas-phase  Fe-DFO-B  (b);  protonated  Fe-DFO-B  in 
solution  using  IEFPCM  (c);  and  protonated  Fe-DFO-B  in  solution  using  72  explicit  FLO  molecules  (d). 
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proven  effective  for  modeling  the  IR  frequencies  of  other 
Fe-organic  complexes  (Domagal-Goldman  and  Kubicki, 
2008),  and  its  use  in  this  study  slightly  reduced  the  disagree¬ 
ment  between  predicted  and  measured  IR  frequencies. 

The  Integral  Equation  Formalism  Polarizable  Contin¬ 
uum  Model  (IEFPCM,  Cances  et  al.,  1997)  was  used  to  ac¬ 
count  for  solvation  implicitly.  When  using  IEFPCM,  the 
large  size  of  the  siderophore  and  corresponding  complexity 
of  the  solvation  cavity  caused  difficulties  with  energy  mini¬ 
mizations.  These  difficulties  were  alleviated  by  using  a  smal¬ 
ler  basis  set  (6-3 1G  instead  of  6-311G(d,p))  and  by 
decreasing  the  number  of  spheres  used  to  create  the 
IEFPCM  cavity.  The  latter  modification  was  accomplished 
by  increasing  the  minimum  radii  of  the  spheres  from  0.2  A 
to  0.5  A  and  reducing  the  default  overlap  index  in  Gaussian 
03  (Frisch  et  al.,  2004)  from  0.89  to  0.8. 

Explicit  solvation  was  also  used  to  model  the  Fe(III)- 
DFO-B  complex  in  a  water  droplet  composed  of  72  H20 
molecules.  After  an  initial  “water  droplet”  structure  was 
obtained  from  a  force-field  based  energy  minimization  in 
Cerius2  (Accelrys,  2003),  a  DFT  energy  minimization  of 
the  complex  was  performed  using  the  CEP-121G  (Walker 
et  al.,  1983;  Cundari  and  Stevens,  1993)  basis  set  for  Fe, 
and  the  3-21G(d,p)  basis  set  for  the  remaining  atoms.  A  lar¬ 
ger  basis  set  for  the  C,  N,  O,  and  H  atoms  was  impractical 
considering  the  system  size  and  availability  of  computa¬ 
tional  resources. 

To  examine  the  contributions  of  covalent  bonding  to  the 
high  chemical  specificity  for  Fe  by  DFO-B,  the  orbital 
occupancies  of  Fe(II)-,  Fe(III)-,  and  Al(III)-DFO-B  were 
calculated  using  NBO  analysis.  Single-point  NBO  calcula¬ 
tions  were  performed  using  the  B3LYP/6-311++G(d,p) 
method  and  the  NBO  3.0  utility  (Foster  and  Weinhold, 
1980;  Reed  and  Weinhold,  1983;  Reed  et  al.,  1985,  1988; 
Carpenter  and  Weinhold,  1988)  included  with  Gaussian 
03.  The  single-point  NBO  calculations  were  performed  on 
the  B3LYP/6-311G(d)  energy -minimized  geometries. 

2.3.  Calculation  of  isotope  fractionation  factors 

Equilibrium  isotope  fractionation  factors  were  calcu¬ 
lated  from  the  predicted  vibrational  frequencies  of  the  Fe- 
DFO-B  complexes  (56Fe  and  54Fe),  according  to  the  equa¬ 


tions  of  Urey  (1947).  First,  we  predicted  the  equilibrium 
constants  for  isotope  exchange  between  a  complex  and  a 
reservoir  of  the  isotopically  active  element.  This  constant 
is  represented  by  /?  and  is  given  by  the  ratio  of  the  partition 
function  with  isotope  A  ( QA )  to  the  partition  function  with 
isotope  B  ( Qb)\ 


P 


Qa 

Qb 


e-hvAi/2kT 
1  _  e~hvAi/kT 


l  _  e-hvBi/kT' 
e~hvBi/2kT 


(i) 


where  the  product  is  taken  over  the  3N-6  vibrational  modes 
(v;),  N  is  the  number  of  atoms  in  the  complex,  k  is  the  Boltz¬ 
mann  constant,  and  T  is  temperature.  In  this  study,  frac¬ 
tionation  was  predicted  for  the  two  most  common 
isotopes  of  Fe,  56Fe  and  54Fe  (A  =  56  and  B  =  54).  This 
form  of  Eq.  (1)  utilizes  the  ratio  of  vibrational  frequencies 
(Redlich,  1935)  to  convert  /?  into  a  reduced  partition  func¬ 
tion.  The  use  of  this  form  of  the  equation  is  critical  when 
using  implicit  solvation  techniques  (Domagal-Goldman 
and  Kubicki,  2008)  and  has  been  shown  to  lead  to  smaller 
error  propagation  from  inaccurate  frequency  predictions 
(Schaad,  1999).  Once  /3  was  calculated,  the  equilibrium  frac¬ 
tionation  factor,  a,  was  calculated  from  the  ratio  of  /3  values 
for  the  two  complexes,  a  (and  by  association,  /?)  are  useful 
in  predicting  equilibrium  fractionations  according  to  the 
following  equation: 

A %Fe(X  -  Y)  «  1000  •  ln(ax_r) 

=  1000  •  [ln(/y  -  ln(/»r)].  (2) 


3.  RESULTS  AND  DISCUSSION 
3.1.  Basis  set  effects:  local  bonding  structure 

In  Table  1,  the  structure  of  the  energy-minimized  gas- 
phase  Fe(III)-DFO-B  complex  is  compared  to  both  experi¬ 
mental  single-crystal  X-ray  diffraction  data  (Dhungana 
et  al.,  2001)  and  solution  EXAFS  measurements  (Edwards 
and  Myneni,  2005).  When  compared  to  the  single-crystal  X- 
ray  diffraction  data,  the  MO/DFT  predictions  overesti¬ 
mated  the  Fe-O(C)  bond  lengths  by  0.06-0.08  A  and 
underestimated  the  Fe-O(N)  bond  lengths  by  about 
0.02  A.  The  average  Fe-O  bond  lengths  compared  well  to 


Table  1 

Fe(III)-DFO-B  structure  based  upon  single-crystal  X-ray  diffraction  measurements  (a — Dhungana  et  al.,  2001),  EXAFS  measurements  of  the 
Fe-DFO-B  complex  in  solution  (b — Edwards  et  al.,  2005),  and  gas-phase  MO/DFT  predictions  (this  study).  The  MO/DFT  predictions  are 
labeled  by  the  basis  set.  Fe-O(C)  and  Fe-O(N)  represent  Fe-O  bond  lengths,  where  the  O  atom  is  bound  to  a  C  or  N  atom.  All  values  are 
averaged  (A  for  bond  lengths  and  °  for  angles). 


Fe(III)-DFO-B 

Parameter 

Exp. 

6-31G 

6-31G(d) 

6-31+G(d) 

6-31  lG(d) 

6-311G(d,p) 

6-311+G(d) 

Fe-O(C) 

2.04a 

2.10 

2.10 

2.11 

2.11 

2.07 

2.12 

Fe-O(N) 

1.98a 

1.98 

1.96 

1.98 

1.96 

1.97 

1.98 

Fe-O  (All) 

2.0  lb 

2.04 

2.03 

2.04 

2.04 

2.02 

2.05 

Fe-C/N  (All) 

2.84b 

2.88 

2.83 

2.84 

2.83 

2.82 

2.85 

C-O 

1.28a 

1.30 

1.27 

1.27 

1.26 

1.28 

1.26 

N-O 

1.38a 

1.41 

1.37 

1.37 

1.36 

1.37 

1.36 

N-C 

1.32a 

1.34 

1.33 

1.33 

1.33 

1.33 

1.33 

O-Fe-O 

78. 7a 

77.6 

77.8 

77.3 

77.3 

77.7 

76.8 

Axial  O-Fe-O 

164.2a 

165.7 

164.9 

164.9 

166.0 

161.8 

164.7 
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the  solution  EXAFS  measurements,  however,  which  is 
important  because  the  average  Fe-O  bond  length  is  a  good 
predictor  of  the  /3  factors  (Domagal-Goldman  and  Kubicki, 
2008).  Each  of  the  other  predicted  bond  lengths  were  within 
0.03  A  of  the  experimentally  measured  values.  Axial  O-Fe- 
O  bond  angles  were  overestimated  by  <2°.  Good  agreement 
between  the  6-31G(d)  and  6-31+G(d)  results  showed  that 
addition  of  diffuse  functions  did  not  significantly  impact 
the  predictions.  Fikewise,  similarity  between  the  6-31G(d) 
and  6-31  lG(d)  results  showed  that  triple-versus  double-zeta 
basis  sets  had  a  modest  effect  on  the  predictions.  In  terms  of 
agreement  with  the  experimental  measurements,  the  most 
cost-effective  basis  set  examined  was  6-31G(d). 

3.2.  Solvation  method  effects:  large-scale  structure 

The  effects  of  solvation  on  the  Fe(III)-DFO-B  structure 
were  also  considered.  Although  solvation  did  not  dramati¬ 
cally  affect  the  local  bonding  environment  surrounding  the 
Fe  atom,  it  did  affect  the  overall  structure  of  the  complex. 
For  example,  when  modeling  the  protonated  (i.e.,  -NH3+) 
gas-phase  complex,  the  amine  group  formed  H-bonds  with 
the  oxygen  atoms  coordinated  to  the  Fe  atom.  Despite  the 
dramatic  effects  of  these  intramolecular  H-bonds  on  the  com¬ 
plex  (Fig.  1),  there  was  not  a  significant  effect  on  the  local 
bonding  environment  surrounding  the  Fe  atom  (Table  2). 

To  model  the  protonated  structure  more  reliably,  we 
used  two  different  methods  to  account  for  solvation.  In 
one  scenario,  hydration  was  accounted  for  explicitly  by  sur- 


Table  2 

Fe(III)-DFO-B  structural  data  from  experimental  measurements 
(Exp.)  and  MO/DFT  predictions.  The  experimental  measurements 
are  from  X-ray  diffraction  measurements  of  crystalline  Fe-DFO-B 
(a — Dhungana  et  al.,  2001),  and  EXAFS  measurements  of  Fe- 
DFO-B  in  solution  (b — Edwards  et  al.,  2005).  The  MO/DFT 
predictions  correspond  to  deprotonated  gas-phase  Fe-DFO-B 
(DFO-B),  protonated  gas-phase  Fe-DFO-B+  (DFO-B+),  proton¬ 
ated  Fe-DFO-B  in  solution  using  the  IEFPCM  method  (PCM), 
and  protonated  Fe-DFO-B  in  solution  using  72  explicit  H20 
molecules  (Droplet).  The  ‘Amine  N-Fe’  distance  is  not  an  actual 
bond  length,  but  represents  the  distance  between  the  N  in  the  amine 
functional  group  of  DFO-B  and  the  Fe  atom.  The  short  amine  N- 
Fe  distance  in  the  DFO-B+  model  (3.74  A)  is  associated  with  intra¬ 
molecular  H-bonding  that  is  an  artifact  of  the  model  and  highlights 
the  need  to  take  into  account  the  effects  of  solvation.  All  values 
except  the  N-Fe  distance  are  averaged  (A  for  bond  lengths  and  ° 
for  angles). 


Fe  (III)-DFO-B 

Parameter 

Exp. 

DFO-B 

DFO-B+ 

PCM 

Droplet 

Fe-O(C) 

2.04a 

2.11 

2.14 

2.08 

2.07 

Fe-O(N) 

1.98a 

1.96 

1.96 

1.99 

1.97 

Fe-O  (All) 

2.01b 

2.04 

2.06 

2.03 

2.02 

Fe-C/N  (All) 

2.84b 

2.83 

2.89 

2.87 

2.88 

C-O 

1.28a 

1.26 

1.31 

1.30 

1.31 

N-O 

1.38a 

1.36 

1.40 

1.41 

1.44 

N-C 

1.32a 

1.33 

1.33 

1.33 

1.32 

O-Fe-O 

78. 7a 

77.3 

76.8 

77.4 

79.2 

Axial  O-Fe-O 

164.2a 

166.0 

163.7 

162.0 

168.1 

Amine  N-Fe 

8.466a 

8.77 

3.74 

8.89 

7.48 

rounding  the  complex  with  72  H20  molecules  (e.g.,  water 
droplet).  In  the  second  scenario,  the  IEFPCM  (Cances 
et  al.,  1997)  method  was  used  to  implicitly  account  for 
hydration  (Table  2).  Differences  in  the  predicted  bond 
lengths  for  these  two  solvation  methods  were  less  than 
0.03  A  for  bonds  directly  involved  in  Fe  complexation. 
Bond  angle  differences  were  1.8°  for  the  O-Fe-O  angles 
and  6.1°  for  the  O-Fe-O  axial  angles.  For  both  of  the  sol¬ 
vated  model  calculations,  the  predicted  Fe-O  bond  lengths 
exhibited  better  agreement  with  experiment  compared  to 
the  gas-phase  calculations  described  in  Section  3.1. 

For  the  solvated  model  calculations,  major  artificial 
changes  to  the  Fe(III)-DFO-B  complex  resulting  from  in- 
tra-molecular  H-bonding  were  absent,  as  evidenced  by  the 
relatively  large  distances  between  the  amine  N  and  Fe 
atoms.  This  distance  changed  by  <1.5  A  between  the  depro¬ 
tonated  gas-phase  model,  the  implicitly  solvated  model,  and 
the  explicitly  solvated  model.  This  difference  is  small  com¬ 
pared  to  the  to  the  A  decrease  in  the  Fe-N  distance  that 
appeared  in  gas-phase  calculations  of  protonated  Fe-DFO- 
B.  Agreement  between  the  deprotonated  gas-phase  model 
and  the  solvated  models  implies  that  the  structure  of  Fe- 
DFO-B  can  be  reasonably  approximated  by  gas-phase 
models  if  artificial  intra-molecular  H-bonding  is  reasonably 
minimized. 

3.3.  Chemical  specificity  of  DFO-B  for  Fe(III) 

In  Table  3,  the  local  bonding  environment  of  uncom- 
plexed  DFO-B  is  compared  with  Fe(III)-,  Fe(II)-,  and 
Al(III)-DFO-B  to  evaluate  differences  that  may  be  qualita¬ 
tively  related  to  the  ligand-metal  affinity.  Replacing  Fe(III) 
with  Fe(II)  resulted  in  lengthening  of  all  the  Fe-O  bonds  by 
more  than  0.1  A  and  a  decrease  in  the  O-Fe-O  angles  by  1- 
2.5°.  Replacing  Fe(III)  with  Al(III)  resulted  in  shortening 
of  the  M-O  bonds  by  0.08-0.15  A  and  an  increase  in  the 
O-M-O  angles  of  3-4°.  The  effects  of  substituting  different 
metals  into  the  DFO-B  complex  were  limited  to  the  local 
bonding  environment.  C-O,  N-O,  and  N-C  bond  length 
changes  were  smaller  than  their  associated  errors,  and  the 
overall  geometry  of  the  complex  was  relatively  unchanged 
(Supplemental  Fig.  1).  Changes  to  the  local  bonding  envi¬ 
ronment  were  significantly  larger  than  differences  arising 
from  use  of  various  basis  sets,  solvation  methods,  and  dis¬ 
agreements  between  predicted  and  experimental  structures. 
The  relative  length  of  the  M-O  bonds  followed  the  order, 
Fe(II)  >  Fe(III)  >  Al(III),  consistent  with  differences  in  the 
ionic  radii  of  the  metal  cations  (0.78  A  for  Fe(II),  0.645  A 
for  Fe(III),  and  0.535  A  for  Al(III)).  However,  the  shorter 
length  of  the  Al(III)-DFO-B  bonds  was  inconsistent  with 
the  preference  of  DFO-B  for  Fe(III)  over  Al(III).  The  abil¬ 
ity  of  DFO-B  to  preferentially  complex  Fe(III)  when 
Fe(III)  and  Al(III)  are  co-present  in  soils  must  therefore 
arise  from  some  other  factor. 

The  differences  between  the  metal-oxygen  bonds  in 
Fe(III)-DFO-B  and  those  in  Fe(II)-DFO-B  and  Al(III)- 
DFO-B  were  revealed  by  NBO  analysis  (Table  4).  NBO 
analysis  assigns  molecular  orbital  electron  densities  using 
a  single  electron  projection  of  the  full  electron  distribution, 
thus  allowing  the  analysis  of  single  electron  orbitals  (Wein- 
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Table  3 


MO/DFT  predictions  of  changes  to  the  local  bonding  environment  of  DFO-B  without  a  metal  and  complexed  with  Fe(II),  Fe(III),  and 
Al(III).  Values  were  averaged  and  standard  deviations  are  listed  in  parentheses  (A  for  bond  lengths  and  °  for  angles).  M-O  bonds  changed 
according  to  the  size  of  the  ionic  radii  of  the  metal  (Fe(II)  >  Fe(III)  >  Al(III)),  but  the  C-O,  N-O,  and  N-C  bonds  were  unchanged. 


Parameter 

Edwardsa 

DFO-B 

Fe(II) 

Fe(III) 

Al(III) 

M-O(C) 

— 

— 

2.24  (0.03) 

2.11  (0.02) 

1.96  (0.01) 

M-O(N) 

2.07  (0.01) 

1.97  (0.01) 

1.89  (0.01) 

C-O 

1.24 

1.23 

1.26  (0.01) 

1.26  (0.00) 

1.27  (0.00) 

N-O 

1.41 

1.41 

1.36  (0.01) 

1.36  (0.00) 

1.37  (0.00) 

N-C 

1.36 

1.36 

1.33  (0.01) 

1.33  (0.00) 

1.32  (0.00) 

O-M-O 

— 

— 

74.8  (0.6) 

77.3  (0.5) 

81.5  (0.5) 

Axial  O-M-O 

— 

— 

165.1  (4.3) 

166.0  (1.2) 

169.8  (1.1) 

a  Edwards  et  al.  (2005). 


hold  and  Landis,  2001).  In  this  study,  NBO  analysis  was 
used  for  quantification  of  the  metal  cation  atomic  charges, 
the  occupancy  of  atomic  orbitals,  and  the  atomic  orbital 
contributions  to  molecular  orbitals.  The  occupancy  of  nat¬ 
ural  atomic  orbitals  agreed  with  the  occupancy  expected  for 
the  given  nuclear  charges  and  calculated  electronic  charges 
of  the  metals.  The  natural  electronic  configurations  deter¬ 
mined  by  the  NBO  analysis  were:  [cor- 

e]4s0'263d5'914p°'014d004  with  a  charge  of  +1.77  for 
Fe(III),  [core]4s°'213d6'264p0'014d0'02  with  a  charge  of 
+1.49  for  Fe(II),  and  [core]3s0363p0'723d°'024p0'01  with  a 
charge  of  +1.88  for  Al(III).  This  qualitatively  agrees  with 
expectations  that  Al(III)  and  Fe(III)  have  greater  positive 
charges  than  Fe(II),  and  with  the  expected  valence  orbital 
occupancies  for  each  metal  (mostly  4s  and  3d  in  Fe  and 
3s  and  3p  in  Al).  The  ~6  d-electrons  for  Fe(III)  is  in  agree¬ 
ment  with  experimental  evidence  for  ligand-to-metal  charge 
transfer  between  DFO-B  and  Fe(III)  (Edwards  et  al.,  2005). 

The  NBO  predictions  indicate  that  Fe(II)-  and  Al(III)- 
DFO-B  bonds  are  completely  ionic,  but  that  four  of  the 
Fe(III)-DFO-B  bonds  are  covalent.  Table  4  shows  the  orbi¬ 
tal  occupancies  for  Fe(II)-,  Fe(III)-  and  Al(III)-DFO-B  for 
the  metal  atoms  in  each  of  these  complexes  and  for  the  four 
O  atoms  that  form  the  covalent  Fe-O  bonds  in  Fe(III)- 
DFO-B,  labeled  in  Fig.  2.  Fe(II)-DFO-B  and  Al(III)- 
DFO-B  did  not  exhibit  any  covalent  bonding  between  the 
metal  and  DFO-B,  as  evidenced  by  the  lack  of  Fe-O  and 
Al-O  participation  by  electrons  from  O  atoms  in  these 
complexes.  This  contrasts  with  the  O  atoms  in  the 
Fe(III)-DFO-B  complex,  for  which  four  covalent  Fe-O 
bonds  were  predicted.  The  presence  of  covalent  bonding  be¬ 
tween  the  sp  orbital  of  O  and  the  sd  orbital  of  Fe(III)  may 
effectively  result  in  the  high  chemical  specificity  of  DFO-B 
for  Fe(III). 

3.4.  IR  frequencies 

The  predicted  IR  frequencies  for  the  de-protonated  gas- 
phase  Fe(III)-DFO-B  complex  and  protonated  Fe(III)- 
DFO-B  in  solution  (modeled  implicitly  using  the  IEFPCM 
method)  reproduced  the  experimentally  measured  IR  fre¬ 
quencies  with  reasonable  accuracy  (Table  5  and  Supple¬ 
mental  Fig.  2,  and  EA-1).  Theoretical  spectra  (Fig.  3) 
were  generated  in  Gaussview  (Dennington  et  al.,  2003)  by 
applying  peaks  with  a  half-width  of  5  cm-1  to  the  scaled 


frequencies.  The  mean  difference  between  the  predicted 
gas-phase  IR  frequencies  presented  here  and  the  experimen¬ 
tal  IR  frequencies  measured  by  Edwards  et  al.  (2005)  was 
19  cm-1  (average  error  of  1.3%,  and  root-mean-squared 
difference  (RMSD)  of  26  cm-1).  The  mean  difference  be¬ 
tween  the  IEFPCM  predicted  IR  frequencies  and  experi¬ 
mental  IR  frequencies  was  15  cm-1  (average  error  of  1.2% 
and  RMSD  of  18  cm-1).  The  accuracy  of  the  IR  frequency 
predictions  was  comparable  to  a  previous  theoretical  study 
of  an  analog  of  DFO-B  (Table  5,  Edwards  et  al.,  2005), 
which  exhibited  an  average  error  23  cm-1  (1.5%  average  er¬ 
ror  and  RMSD  of  32  cm-1). 

The  advantage  of  including  solvation  was  more  appar¬ 
ent  when  examining  the  best  fit  line  through  the  model  pre¬ 
dictions.  If  the  model  predictions  exactly  reproduce  the 
experimental  measurements,  then  all  model  frequencies 
would  fall  on  the  line  y  =  x  when  plotted  as  a  function  of 
experimental  frequencies  (Supplemental  Fig.  2).  The  best 
fit  line  for  the  gas-phase  calculation  was 

y  =  1.0287v  —  38.97,  with  an  r2  value  of  0.98.  The  best  fit 
line  for  the  IEFPCM  model  calculation  was 

y  =  1.001 9x  —  1.513,  with  an  r2  value  0.99.  The  decrease 
in  error  associated  with  the  IEFPCM  method  resulted  pri¬ 
marily  from  better  prediction  of  the  C=0  stretching  modes. 
The  IEFPCM  calculation  of  the  1577  cm-1  and  1627  cm-1 
vibrational  modes  were  in  error  by  5  cm-1  and  28  cm-1, 
respectively,  compared  to  errors  of  34  cm-1  and  61  cm-1 
for  the  gas-phase  calculations,  and  44  cm-1  and  78  cm-1, 
respectively,  for  calculations  of  a  gas-phase  DFO-B  analog 
(Edwards  et  al.,  2005).  The  accuracy  of  the  models  used  in 
this  study  was  an  improvement  over  the  predictions  from 
aHA  analogs  of  DFO-B.  This  is  important  for  calculating 
equilibrium  isotope  constants  as  they  depend  exponentially 
on  the  predicted  frequencies.  However,  the  modest  changes 
in  the  errors  justify  the  methods  employed  by  Edwards 
et  al.  (2005),  as  their  methods  were  much  less  computation¬ 
ally  demanding.  Nevertheless,  their  study  did  not  include 
Fe  isotope  fractionations. 

The  root-mean-squared  error  in  the  predicted  IR  fre¬ 
quencies  from  this  study  was  26  cm-1,  and  the  root-mean- 
squared  error  in  the  predicted  frequencies  for  various  Fe- 
organic  complexes,  including  those  from  a  previous  study 
(Domagal-Goldman  and  Kubicki,  2008),  was  31  cm-1.  Be¬ 
cause  of  the  lack  of  low-frequency,  isotopically-sensitive 
vibrational  modes  in  the  DFO-B  measurements  (Edwards 
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Table  4 

NBO  analysis  of  DFO-B  complexes  with  Fe(II),  Fe(III),  and  Al(III).  Single-point  calculations  were  performed  with  the  B3LYP/6-311++G(d,p)  method.  Orbital  occupancies  are  shown  for  the 
metal  atoms  and  the  4  O  atoms  that  exhibit  covalent  bonding  with  Fe(III).  By  comparing  the  orbital  occupancies  of  similar  atoms  in  different  complexes,  the  source  of  the  electrons  participating  in 
covalent  bonds  can  be  tentatively  identified.  Integers  next  to  the  O  atoms  are  strictly  for  labeling  purposes  (see  Fig.  2).  2-center  covalent  bonds  are  labeled  A-B,  where  A  and  B  are  the  two  atoms 
donating  electrons  to  the  bond.  A-B*  indicates  an  anti-bonding  orbital  between  A  and  B.  Lone  pair  valence  electrons  are  labeled  LP.  LP*  indicates  an  unoccupied  lone  pair  orbital.  Contributions 
from  individual  atomic  orbitals  are  also  listed.  Note  the  only  complex  for  which  the  metal  cation  has  electrons  in  covalent  bonds  is  Fe(III)-DFO-B. 


Atom 

Fe(II) 

Fe(III) 

Al(III) 

Type 

Occ. 

Composition 

Type 

Occ. 

Composition 

Type 

Occ. 

Composition 

O  27 

O-N 

0.9931 

9%  Sq,  40%  po,  10%  Sn,  30%  pN 

O-N 

0.9917 

9%  Sq,  42%  po,  12%  Sn,  37%  pN 

O-N 

0.9938 

9%  s0,  42%  po, 

12%  sn,  37%  pN 

LP 

0.9842 

59%  s,  42%  p 

O-Fe 

0.9813 

17%  So, 70%  po,4%  Spe,  9%  dpe 

LP 

0.9821 

36%  s,  64%  p 

LP 

0.9211 

21%  s,  79%  p 

O-Fe 

0.9530 

89%  Po,  n%  dFe 

LP 

0.9397 

100%  p 

LP 

0.9128 

3%  s,  97%  p 

LP 

0.9794 

62%  s,  38%  p 

LP 

0.9094 

46%  s,  54%  p 

O  28 

O-N 

0.9931 

9%  Sq5  40%  po,  13%  Sn,  38%  pN 

O-N 

0.9918 

9%  Sq,  42%  po,  12%  Sn,  37%  pN 

O-N 

0.9938 

9%  s0,  42%  po, 

12%  Sn,  37  pN 

LP 

0.9832 

61%  s,  39%  p 

O-N 

0.9429 

65%  po,  35%  pN 

LP 

0.9818 

36%  s,  64%  p 

LP 

0.9208 

7%  s,  93%  p 

O-N* 

0.6967 

35%  po,  65%  pN 

LP 

0.9414 

100%  p 

LP 

0.9134 

15%  s,  85%  p 

O-Fe 

0.9817 

17%  Sq,  70%  po,  4%  Spe,  9%  dpe 

LP 

0.9100 

46%  s,  54%  p 

LP 

0.9794 

63%  s,  37%  p 

O  57 

O-N 

0.9750 

52%  po,  48%  pN 

O-C 

0.9934 

24%  s0,  41%  po,  11%  sc,  24%  pc 

O-C 

0.9939 

25%  s0,  40%  po,  1 1%  sc,  24%  pc 

O-N 

0.6130 

48%  po,  52%  Pn 

O-Fe 

0.9640 

12%  Sq,  79%  po,  3%  Spe,  6%  dpe 

LP 

0.9704 

36%  s,  64%  p 

O-C 

0.9947 

25%  s0,  39%  po,  11%  sc,  25%  pc 

LP 

0.9744 

50%  s,  50%  p 

LP 

0.909 

26%  s,  74%  p 

LP 

0.9768 

55%  s,  45%  p 

LP 

0.7985 

100%  p 

LP 

0.8352 

100%  p 

LP 

0.9184 

6%  s,  94%  p 

O  72 

O-N 

0.9935 

9%  Sq,  40%  po,  13%  S]\r,  38%  pN 

O-Fe 

0.9512 

87%  po,  13%  sd25'64^ 

O-N 

0.9938 

9%  s0,  42%  po, 

12%  Sn,  37%  pN 

LP 

0.9839 

60%  s,  40%  p 

O-N 

0.9915 

10%  Sq,  41%  po,  12%  Sn,  37%  pN 

LP 

0.9819 

36%  s,  64%  p 

LP 

0.9245 

13%  s,  87%  p 

LP 

0.9786 

62%  s,  38%  p 

LP 

0.9399 

100%  p 

LP 

0.9076 

9%  s,  91%  p 

LP 

0.8623 

18%  s,  82%  p 

LP 

0.9087 

46%  s,  54%  p 

Fe/Al 

* 

LP 

0.1241 

100%  d 

Fe-027 

0.9813 

17%  Sq,  70%  po,  4%  Spe,  9%  dpe 

* 

LP 

0.1780 

100%  s 

LP* 

0.1209 

100%  d 

Fe-027 

0.9813 

89%  Po,  n%  dFe 

LP* 

0.1232 

100%  p 

LP* 

0.1011 

98%  s,  2%  d 

Fe-028 

0.9817 

17%  Sq,  70%  po,  4%  Spe,  9%  dpe 

LP* 

0.1225 

100%  p 

LP* 

0.0658 

100%  d 

Fe-057 

0.9640 

12%  Sq,  79%  po,  3%  Spe,  6%  dpe 

LP* 

0.1166 

100%  p 

LP* 

0.0569 

100%  d 

Fe-072 

0.9512 

87%  po,  13%  dFe 

LP 

0.9726 

100%  d 

LP 

0.1577 

100%  d 
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Fig.  2.  Close  up  of  the  Fe-O  bonds  from  the  gas-phase,  depro- 
tonated  Fe-DFO-B  optimized  geometry  (Fig.  lb).  O  atoms  (in  red) 
bound  to  the  Fe  atom  are  labeled  with  numbers  consistent  with 
those  listed  in  Table  2.  The  numbering  is  consistent  across  different 
models — Fe(II),  Fe(III),  and  Al(III).  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  paper). 

et  al.,  2005),  we  use  the  latter  value  when  evaluating  the 
accuracy  of  the  Fe  isotope  fractionation  factors.  The  errors 
of  the  predicted  IR  frequencies  were  translated  into  errors 
in  predictions  of  equilibrium  isotope  fractionation  factors 
by  randomly  generating  20  errors  for  each  predicted  fre¬ 
quency,  using  a  standard  deviation  of  31  cm-1.  This  al¬ 
lowed  the  /3  factors  to  be  re-calculated  20  times  with 


Fig.  3.  Predicted  IR  spectra  for  B3LYP/6-311G(d)  gas-phase 
calculations  of  Fe(III)-DFO-B  (vacuum)  and  B3LYP/6-31G(d) 
IEFPCM  calculations  (IEFPCM).  Vibrational  frequencies  were 
calculated  as  described  in  Section  2.2,  and  the  spectra  were 
generated  from  these  frequencies  using  a  full-width  half-maximum 
of  5  cm-1  in  Gaussview  (Dennington  et  al.,  2003). 

random  errors  applied  to  each  vibrational  mode.  The  con¬ 
fidence  intervals  reported  in  this  study  represent  the  stan¬ 
dard  deviations  of  the  20  re-calculations  of  /3 .  Errors  in 
predictions  of  1000-ln(a)  were  calculated  by  taking  the 
square  root  of  the  sum  of  the  squares  of  the  1000-ln(^)  er¬ 
rors,  and  were  ~0.25%o.  A  more  detailed  discussion  on  this 
and  other  sources  of  error  in  MO/DFT  predictions  of  Fe 
isotope  fractionation  factors  can  be  found  in  previous  work 
(Domagal-Goldman  and  Kubicki,  2008). 

3.5.  Isotopic  equilibrium  in  solution 

Accurate  prediction  of  the  structures  and  IR  vibrational 
frequencies  for  Fe-DFO-B  complexes  enabled  the  calcula¬ 
tion  of  equilibrium  isotope  fractionations.  1000-ln ([3)  values 


Table  5 

Comparison  of  observed  and  calculated  IR  vibrational  frequencies  (scaled).  Frequencies  are  listed  in  cm-1.  Intensities  are  listed  in  cm_2/mol/ 
L,  converted  from  km/mol  using  a  factor  of  100  (Kubicki  et  al.,  1993).  The  calculations  from  Edwards  et  al.  (2005)  were  for  an  uncomplexed 
model,  and  those  from  this  study  are  for  the  Fe(III)-DFO-B  complex.  The  assignments  listed  are  from  IR  measurements  published  by 
Edwards  et  al.  (2005).  All  predicted  frequencies  from  this  study  (including  those  that  do  not  correlate  to  experimentally-determined 
frequencies)  are  listed  in  Electronic  Annex  EA-1,  along  with  their  intensities  and  56v/54v  ratios. 


Edwards  et  al.  (2005) 

This  study 

Exp. 

Theoretical 

Assignment 

Gas-Phase 

Inten. 

IEFPCM 

Inten. 

1049 

1036 

vC-C 

vC-N 

1046 

2205 

1035 

2010 

1184 

1188 

6C-H  twist 

5C-H  wag 

1196 

4045 

1170 

7672 

1209 

— 

— 

— 

— 

1236 

6425 

1355 

1351 

5C-H  wag 

6C-H  twist 

1324 

5078 

1347 

3933 

1374 

1375 

vC-N 

6C-H  wag 

5C-H  rock 

6C-H  twist 

1380 

2806 

1390 

20273 

1421 

1427 

5CCN  bend 

5CCC  bend 

1418 

4221 

1428 

4969 

1453 

1486 

5NOH  bend 
vC-N  rock 
vC-N  rock 

1449 

5642 

1477 

7276 

1577 

1621 

vC=0 

1543 

23114 

1582 

32335 

1627 

1705 

vC=0 

vC-N 

1688 

45111 

1599 

10789 
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have  been  predicted  for  both  Fe(II)  and  Fe(III)  coordinated 
with  H20  (Jarzecki  et  al,  2004;  Anbar  et  al.,  2005),  com- 
plexed  by  oxalate  and  catecholate  (Domagal-Goldman 
and  Kubicki,  2008),  and  now  complexed  by  DFO-B  (Table 
6  and  Fig.  4).  Equilibrium  Fe  isotope  fractionations  can 
thus  be  predicted  for  these  different  Fe  species  in  solution 
(see  Eq.  (2))  by  taking  the  differences  of  their  1000-ln (/?)  val¬ 
ues  (i.e.,  by  calculating  the  1000-ln(a)  values  for  the  reaction 
in  question). 

We  predict  an  equilibrium  fractionation  of 
—0.34  ±  0.25%o  between  Fe(III)-DFO-B  and  Fe(III)- 
(H20)6  at  22  °C  (Table  6),  in  disagreement  with  the 
+0.60  dz  0.1 5%0  value  measured  by  Dideriksen  et  al. 
(2008)  between  Fe(III)-DFO-B  and  inorganic  phases  (pri¬ 
marily  composed  of  Fe(III)-(H20)6).  If  the  experimental 
measurements  are  accurate  within  reported  uncertainties, 
either  the  model  predictions  of  the  lOOO  ln(^)  values  for 
the  organic  fraction  are  too  large  or  the  model  predictions 
of  the  1000-ln (/?)  values  for  the  inorganic  fraction  are  too 
small  (or  both,  with  some  combination  of  the  two  sources 
of  error). 

Previous  predictions  of  1000-ln(a)  values  for  equilibrium 
between  inorganic  solution  and  mineral  phases  (Domagal- 
Goldman  and  Kubicki,  2008)  also  differ  from  experimental 
measurements  by  ~1.0%o,  with  the  predictions  of  aqueous 
phase  lOOO-ln(^)  values  too  large  compared  to  mineral 
phase  values.  Thus,  the  discrepancies  with  experiment  in 
both  theoretical  studies  could  be  caused  by  overestimates 
of  aqueous  phase  inorganic  1000-ln(^).  This  could  be  due, 
in  part,  to  not  accounting  for  all  inorganic  phases  present 
in  an  experiment.  To  test  this  possibility,  we  calculated 
the  1000-ln  (/?)  values  of  Fe(III)-(H20)5(0H),  a  major  com¬ 
ponent  of  inorganic  Fe(III)  in  experiments,  and  of  Fe(III)2- 
(H2O)i0(OH)2,  an  Fe-dimer  complex  (Fig.  5)  that  is  stoi- 
chiometrically  equivalent  to  Fe(III)-(H20)5(0H)  and  may 
be  the  structure  that  exists  in  solution.  The  predicted 
1000-ln(a)  value  for  the  Fe-dimer  complex  at  22  °C  was 
8.20  +  0.12.  This  value  leads  to  a  prediction  of 
+0.45  ±  0.23  for  isotopic  equilibrium  between  Fe(III)- 
DFO-B  and  Fe(III)2-(H20)6(0H)2,  within  experimental  er- 


9 


Fig.  4.  1000-ln(/i)  values  plotted  against  1/T2  (106K-2).  56Fe/54Fe 
increases  towards  the  top  of  the  figure,  and  temperature  increases 
towards  the  left  side  of  the  figure.  Note  the  separation  between 
Fe3+  and  Fe2+  phases,  which  indicates  the  importance  of  redox 
reactions  in  isotopic  fractionations.  Values  for  Fe-(H20)6  and  Fe- 
(oxalate)3  are  from  Domagal-Goldman  and  Kubicki  (2008). 


ror  of  the  isotopic  equilibrium  between  Fe(III)-DFO-B  and 
inorganic  precipitates.  This  agreement  may  suggest  that 
fractionation  is  controlled  by  isotopic  equilibrium  between 
Fe(III)-DFO-B  and  Fe(III)2-(H2O)i0(OH)2  prior  to  precip¬ 
itation  of  inorganic  Fe(III). 

Given  the  degree  to  which  equilibrium  between  Fe(III)- 
DFO-B  and  all  inorganic  phases  was  obtained  in  the  exper¬ 
iments,  and  the  short  timescales  of  precipitation,  it  is  unli¬ 
kely  that  equilibrium  between  Fe(III)-DFO-B  and  a  single 
inorganic  phase  is  responsible  for  the  measured  fraction¬ 
ation.  A  better  approach  would  be  to  weigh  the  lOOO-ln(^) 
values  of  the  inorganic  phases  by  their  reported  abundances 


Table  6 

Predicted  1000-ln(^)  values  for  Fe(III)  complexed  with  water,  oxalate,  catecholate,  hydroxide,  and  DFO-B.  These  calculations  used  the 
UB3LYP/6-311G(d,p)  methodology,  IEFPCM  to  handle  solvation  effects,  and  a  frequency  scale  factor  of  0.956. 


7TC) 

Fem(H20)63+ 

Fenl(0H)(H20)52+ 

(Fem)2(0H)2(H20)62+ 

FemDFOB+ 

Fem(Ox)33- 

Fenl(Cat)33 

0 

9.98  +  0.18 

10.19  +  0.18 

9.48  +  0.14 

9.62  ±  0.22 

9.09  +  0.18 

8.84  +  0.25 

10 

9.34  +  0.17 

9.54  +  0.17 

8.86  +  0.13 

8.99  +  0.21 

8.50  +  0.17 

8.26  ±  0.23 

20 

8.76  +  0.16 

8.96  +  0.16 

8.31+0.12 

8.42  ±  0.20 

7.96  +  0.16 

7.73  ±  0.22 

22 

8.65  +  0.16 

8.84  +  0.16 

8.20  +  0.12 

8.31+0.20 

7.86  +  0.16 

7.63  ±  0.22 

25 

8.49  +  0.16 

8.68  +  0.15 

8.05  +  0.12 

8.15  +  0.19 

7.71+0.16 

7.49  ±  0.21 

30 

8.23  +  0.15 

8.42  +  0.15 

7.80  +  0.12 

7.90  +  0.18 

7.47  +  0.15 

7.26  +  0.21 

40 

7.74  +  0.15 

7.93  +  0.14 

7.34  +  0.11 

7.43  +  0.17 

7.02  +  0.14 

6.82  +  0.19 

50 

7.30  +  0.14 

7.48  +  0.13 

6.92  +  0.11 

7.00  +  0.16 

6.61+0.13 

6.42  +  0.18 

60 

6.89  +  0.13 

7.07  +  0.13 

6.53  +  0.10 

6.60  +  0.15 

6.24  +  0.13 

6.06  +  0.17 

70 

6.52  +  0.13 

6.69  +  0.12 

6.17  +  0.09 

6.24  +  0.15 

5.90  +  0.12 

5.72  +  0.16 

80 

6.17  +  0.12 

6.34  +  0.11 

5.84  ±  0.09 

5.90  +  0.14 

5.58  +  0.11 

5.41+0.16 

90 

5.85  +  0.11 

6.02  +  0.11 

5.54  ±  0.09 

5.59  +  0.13 

5.29  +  0.11 

5.13  +  0.15 

100 

5.56  +  0.11 

5.72  +  0.10 

5.26  ±  0.08 

5.31+0.13 

5.02  +  0.10 

4.87  +  0.14 

200 

3.52  ±  0.07 

3.63  ±  0.07 

3.33  +  0.05 

3.35  +  0.08 

3.17  +  0.07 

3.07  ±  0.09 

300 

2.42  ±  0.05 

2.51+0.05 

2.29  ±  0.05 

2.30  ±  0.06 

2.18  +  0.05 

2.11  ±0.09 
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Fig.  5.  Structure  of  Fe(III)2-(OH)2-(H20)5,  as  predicted  by 
B3LYP/6-311G(d,p)  utilizing  IEFPCM  to  model  solvation  effects. 
Orange  atoms  represent  Fe(III),  red  atoms  represent  O,  and  white 
atoms  represent  H.  This  is  a  potential  structure  for  Fe(III)  in 
solution,  and  may  help  describe  the  measured  equilibrium  isotope 
fractionation  between  inorganically  bound  Fe  and  Fe  bound  to 
organic  ligands.  (For  interpretation  of  the  references  to  colour  in 
this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
paper.) 

based  upon  the  results  of  thermodynamic  modeling.  The 
abundances  in  the  experiment  were  70%  Fe(III)-(H20)6 
(1000-ln(j8)  =  8.65  +  0.16),  29%  Fe(III)2-(H20)8(0H)2, 
(1000-ln(jff)  =  8.20  ±0.12),  and  1%  Fe(II)-(H20)6, 
(1000-ln(/?)  =  5.66  ±  0.16).  Thus,  the  predicted  1000-ln(jff) 
value  for  the  inorganic  phase  is 
0.70  x  8.65  +  0.29  x  8.20  +  0.01  x  5.66  =  8.46,  leading  to 
a  Fe(III)-DFO-B-Fe-inorganic  fractionation  prediction  of 
—0.15  ±  0.24%o.  This  represents  an  improvement  over  the 
predictions  based  solely  on  Fe(III)-(H20)6,  but  is  not  within 
error  of  the  experimental  measurement,  indicating  that 
either  isotopic  equilibrium  in  experiments  is  dominated  by 
isotopic  exchange  between  Fe(III)-DFO-B  and  Fe(III)2- 
(H2O)i0(OH)2  prior  to  precipitation  of  Fe-oxyhydroxides, 
or  there  is  some  other  source  of  error  in  predicted  or  mea¬ 
sured  1000-ln(a)  values. 

The  predicted  equilibrium  fractionations  between 
Fe(III)-DFO-B  and  Fe(III)-(oxalate)3  (+0.45  ±  0.25%o,  Ta¬ 
ble  6)  and  between  Fe(III)-DFO-B  and  Fe(III)-(catecho- 
late)3  (+0.68  +  0.25%o,  Table  6)  were  close  to  differences 
between  the  isotopic  composition  of  Fe  in  cells  and  Fe 
remaining  in  solution  during  hornblende  dissolution,  mea¬ 
sured  to  be  +0.85  and  +1.17  +  0.20%o  (Brantley  et  al., 
2001).  This  agreement  is  consistent  with  isotopic  equilib¬ 
rium  between  Fe  complexed  with  small  organic  ligands 
and  Fe  complexed  with  DFO-B  and  no  isotopic  fraction¬ 
ation  between  Fe(III)-DFO-B  and  cell  material.  While 
other  fractionations  may  exist  in  such  a  system,  including 
reactions  inside  cells  and  during  transfer  of  Fe  from  sidero- 
phores  to  the  cell,  the  fractionation  reported  by  Brantley 
et  al.  (2001)  can  be  explained  by  the  predicted  equilibrium 
fractionations  between  Fe  complexed  to  small  organic  li¬ 
gands  and  Fe  complexed  to  siderophores. 

Experiments  have  shown  that  the  rate  at  which  Fe  is  lea¬ 
ched  from  mineral  surfaces  increases  when  DFO-B  is  co¬ 


present  with  oxalate  (Cheah  et  al.,  2003).  This  implies  that 
Fe  sequestration  occurs  through  a  multi-step  process:  Fe  is 
first  leached  by  small  organic  ligands,  complexed  by  sidero¬ 
phores  in  solution,  and  delivered  to  cells  via  Fe-sidero- 
phore  complexes.  Our  results  support  this  conceptual 
model  because  the  predicted  equilibrium  Fe  isotope  frac¬ 
tionation  between  Fe(III)-DFO-B  and  Fe(III)-(oxalate)3 
matches  the  experimental  differences  between  the  isotopic 
compositions  of  cells  and  solution  during  hornblende  disso¬ 
lution  (Brantley  et  al.,  2001).  Such  agreement  is  consistent 
with  delivery  of  Fe  to  cells  exclusively  via  Fe(III)-DFO-B 
complexes.  Fractionations  have  also  been  shown  to  exist 
during  dissolution  of  minerals  (Brantley  et  al.,  2001,  2004; 
Wiederhold  et  al.,  2006)  and  during  precipitation  of  Fe(III) 
phases  (Welch  et  al.,  2003).  Thus,  there  may  be  at  least 
three  sources  of  fractionation  during  biological  uptake  of 
Fe:  dissolution  of  the  initial  mineral  phase  by  small  organic 
acids,  equilibrium  between  various  Fe-ligand  complexes  in 
solution,  and  re-precipitation  of  inorganic  mineral  phases. 

4.  CONCLUSIONS 

The  prediction  of  isotope  fractionation  factors,  placed  in 
the  context  of  experimental  measurements,  further  develops 
a  molecular-scale  understanding  of  the  transport  of  Fe 
from  minerals  to  bacterial  cells.  Small  organic  acids  remove 
Fe  from  the  mineral  surface,  and  then  Fe  is  complexed  by 
siderophores  that  ultimately  transport  the  Fe  to  bacterial 
cells.  This  conceptual  model  is  supported  by  our  explana¬ 
tion  of  previous  measurements  of  a  fractionation  between 
Fe  in  solution  and  Fe  in  bacterial  cells  during  the  leaching 
of  Fe  from  hornblende  (Brantley  et  al.,  2001).  The  explana¬ 
tion  of  these  data  is  based  on  our  predicted  fractionation 
factor  for  isotopic  equilibrium  between  Fe(III)-DFO-B 
and  either  Fe(III)-(catecholate)3  or  Fe(III)-(oxalate)3. 

The  role  of  siderophores  in  this  process  is  to  prevent  Fe 
precipitation  by  efficiently  sequestering  it  in  Fe-organic 
complexes.  This  is  only  possible  because  of  the  chemical 
specificity  siderophores  have  for  Fe(III),  which  NBO  anal¬ 
ysis  suggests  is  the  result  of  covalent  metal-oxygen  bonds 
that  are  absent  in  Al(III)-DFO-B  and  Fe(II)-DFO-B.  In 
Fe(III)-DFO-B,  the  Fe(III)-0  bonding  orbitals  are  fully 
populated.  These  orbitals  allow  siderophores  to  efficiently 
sequester  Fe  and  not  other  metals,  thereby  giving  sidero¬ 
phores  a  unique  and  important  role  in  biological  Fe 
acquisition. 
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